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Using the factorization approach of quantum mechanics, we obtain a family of 
isospectral scalar potentials for power law inflationary cosmology. The construction 
is based on a scattering Wheeler-DeWitt solution. These iso-potentials have new 
features, they give a mechanism to end inflation, as well as the possibility to have new 
inflationary epochs. The procedure can be extended to other cosmological models. 
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I. INTRODUCTION 

One of the most active areas of research nowadays is inflationary cosmology, this theoret- 
ical framework solves many classical problems of Standard Big Bang Cosmology. Recently, 
observations have confirmed its predictions of a fiat universe with a nearly scale invariant 
perturbations spectrum. The idea is to introduce a scalar field (spin-0 boson) and a scalar 
potential V{(f)) which encodes in itself the (non-gravitational) self-interactions among the 
scalar particles. This type of models, have also been used within the so called canoni- 
cal Quantum Cosmology (QC) formalism, which deals with the early epoch of the cosmos. 
Scalar fields act as matter sources, and then play an important role in determining the 
evolution of the early universe, where the quantum fluctuations are the seeds for structure 
formation. 

Moreover, these models have appeared in String Theory, in particular in connection to 
the so called string theory landscape as well as in the study of tachyon dynamics. For 

*Electromc address: |alberto@fisica.ugto.inx| 
'''Electronic address: |wguzman@flsica.ugto.nix| 
•l- Electronic address: |nisabido@fisica. ugto.mx| 
^Electronic address: |socorro@fisica.ugto.mx| 



2 



the String Theory landscape (Kobakhidze and Mersini-Houghton, 2004; Douglas, 2003; 
Susskind, 2003), the scalar potential V(0) is usually thought as having many valleys, which 
represent the different vacua solutions, the hope is that the statistics of these vacua could 
explain, for example, the smallness of the cosmological constant (the simplest candidate 
for dark energy). For tachyon dynamics in the unstable D-brane scenario, the scalar po- 
tential for the tachyon effective action around the minimum of the potential has the form 
V(0) = e-"'^/^ (Sen, 2002,2003). This leads to the study of tachyon driven cosmology (Gorini 
2004; Garcia-Compean, Garcia- Jimenez, Obregon and Ramirez, 2005). 

On another front, the study of eigenvalue problems associated with second-order differ- 
ential operators found a renewed interest with the application of the factorization technique 
and its generalizations (Cooper, Khare and Sukhatme, 1995; Fernandez, 1984; Mielnik, 1984; 
Nieto, 1984; Gelfand and Levitan, 1955). SUSY-QM may be considered an equivalent for- 
mulation of the Darboux transformation method, which is well-known in mathematics from 
the original paper of Darboux (Darboux, 1982; Ince, 1926). An essential ingredient is a 
differential operator (Bagrov and Samsonov, 1995) which intertwines two hamiltonians and 
relates their eigenfunctions. When this approach is applied in quantum theory it allows to 
generate a family of exactly solvable local potential starting with a given exactly solvable one 
(Cooper, Khare and Sukhatme, 1995). In nonrelativistic one-dimensional supersymmetric 
quantum mechanics, the factorization technique was applied to the q = factor ordered 
WDW equation corresponding to the FRW cosmological models without matter field ( Rosu 
and Socorro, 1998), where a one-parameter class of strictly isospectral cosmological FRW so- 
lutions was exactly found, representing the wave functions of the universe for that case, also 
iso-spectral solutions for a one-parameter family of closed, radiation-filled FRW quantum 
universe, and for a perfect fluid with barotropic state equation and cosmological constant 
term, for any factor ordering were found (Rosu and Socorro, 1996; Socorro, Reyes and Gel- 
bert, 2003). In this formalism, the family of iso-potentials and wave functions are build with 
respect to a parameter 7j for which we choose the domain [0,oo]. The shape of the wave 
function in the corresponding coordinates is obtained via the "evolution" of the iso-wave 
function when this parameter tends to oo. 

The main purpose of this paper is to apply the Darboux transformation method to obtain 
a family of iso-potentials, to the potential 1^(0) = e""*^/^, which appears in inflationary 
cosmology. To reach this goal, we shall make use of the strictly isospectral scheme based 
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on the general Riccati solution (Cooper, Khare and Sukhatme, 1995; Fernandez, 1984; 
Mielnik, 1984; Nieto, 1984), which is also known as the double Darboux method. This 
scheme has been applied from classical and quantum physics (Mielnik, 1984) to relativistic 
models (Samsonov and Suzko, 2003). This technique has been known for a decade in one- 
dimensional super symmetric quantum mechanics (SUSY-QM) and usually requires nodeless, 
normalizable states of a Schrodinger-like equation. However, Pappademos, Sukhatme, and 
Pagnamenta (1993) showed that the strictly isospectral construction can also be performed 
on non-normalizable states. The resulting potentials have interesting features, in particular 
they solve the one problem the exponential potential has, the lack of a mechanism to end 
inflation. 

This work is organized as follows. In section II we present the classical action with 
the corresponding contributions, this action includes a gravitational part Sg, and S^ for the 
scalar field; also we present the standard quantum scheme with the quantum solution, which 
plays an important role in the isospectral solutions. In section III we review the factorization 
approach of (Susy-QM) and apply it to the inflationary model. Finally section IV is devoted 
to conclusions and outlook. 



II. THE STANDARD QUANTUM SCHEME 

We start with the line element for a homogeneous and isotropic universe, the so called 
Friedmann-Robertson- Walker (FRW) metric, in the form 

dr2 



ds' = -N'(t)dt' + e 



2a(t) 



1 -kr2 



+ r^(dr + sin'M(^' 



(1) 



where a(t) = e"*^*^ is the scale factor, N(t) is the lapse function, and k is the curvature 
constant that takes the values 0, +1, —1, which correspond to a flat, closed or open universe, 
respectively. 

The effective action we will be working, is ( W. Guzman, Sabido, Socorro and Arturo 
Urena Lopez, 2005) 



Stot — Sg + S^ 



(2) 



is a scalar field endowed with a scalar potential V(0) = Vqc ^ 
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The Lagrangian for a FRW cosmological model is 



,3q 



(3) 



At this point, we consider a flat universe (k = 0) 

>.2 X ^2 



(4) 



The canonical momenta are found to be 

« = Y^e-3"n„ , (5a) 
= -Ne-^"n^ . (5b) 

We are now in position to write the corresponding canonical Hamiltonian (Ryan, 1972) 



n = ^ = 12e3"^ 
" da N 

n. = — = -e^"^ 
86 N' 



-^canonical = H^d + U^cf) — NH, 

where Ti is the classical Hamiltonian function, having the following structure 

n = — e-3" [n^ - 12Ul - 24e^"V 

OA L " <P 



(6) 



(7) 



and performing the variation of (jH)) with respect to N, dC/dN = 0, implies the well-known 
result Ti. = 0. The Wheeler-DeWitt (WDW) equation for this model is achieved by replacing 
HqM by — iSqM in Eq. ((Zj); here q'^ = {a,(f)). 

Under a particular factor ordering the WDW reads 



-3a 



H 



24 



da' 



rj2 

12— - 24e^°V(0) 



^ = 0. 



(8) 



or 



□ ^-24e''"V(0)^ = O, (9) 

with d) = -t=, is called the wave function of the universe, □ = —d^ + dj is the two 
dimensional d'Alambertian operator in the q^^ coordinates. From now on we flx the poten- 
tial to V{(f)) = e~'^'^. Applying the factorization method in these variables is technically 
cumbersome, this can be simplifled if we make the following coordinates transformation, 

6 ~ 



X = 6« — A0, y = 



A 



(10) 
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the WDW equation Q takes the form 

d^_^d^_J4V^ (11) 

and by separation variables, \E' = X(x)Y(y) with y = Ay, we obtain the set of differential 
equation for the functions X and Y 

d^X / 71^\ 

+ -/5e" + -L X = 0, 



dx2 V 4 



where 



d^Y 71^ 

d^ + - °' 



3 - (13) 

^"A2-36' ^ 
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we choose for simplicity ^ as a separation constant. The solutions for these equations are 



X(x) = Z±i, (±2i^e''/^j , (14) 
Y(y) = Aoe^^^ + Aie-'^^, (15) 

with Z-ti^ are generic Bessel Functions with pure imaginary order, the wave function is 

^,(x, y) = e^'^^^Zi, (±2i v^e^/2) . (16) 

Since these solutions have the dependence in the parameter rj, the general solution can 
be put as 

^gen = [ G{r])^^dr], (17) 



where G{ri) represents a weigh function. 

The selection of the value of A in (fT^. gives the structure for the Z-ti^, for A > 6, we have 
the modified Bessel function, and for < A < 6, Z-ti^ become the ordinary Bessel function. 
We are now in a position to form a normalizable Gaussian state as a superposition of the 
eigenfunctions (|T6|l . With this in mind a wave packet can be constructed (Kiefer, 1988, 1990). 

We can have different solutions that depend on the value of the parameter A. For A > 6, 
the wave packet can be constructed using the modified Bessel function (see Gradshteyn and 
Ryzhik, 1980 ) 

v&(x, y) = cos (^^y^ K;, (±2v^e^/2j dr] = ^Exp -2 v^e^/^^osh Qy^ . (18) 
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In the range < A < 6 ( which includes the inflation), we can also construct a wave packet, 
but this time using the ordinary Bessel functions 



oo g7rx/2 



sinh(7rx) 



cos ( ) Ji, (±2v^e^/2) dr] = -iExp 



i2V^e"/2cosh ( ^y 



(19) 



We have been using the new variables x and y. Let us now extract some information from 
the semiclassical behavior as a check of our quantum model. The classical solutions can be 
obtained using the semiclassical analysis (WKB-like method). For this, one considers the 
ansatz on the wave function 



and the conditions 



dsy 

da J 



da? 



from this, the Einstein-Hamilton- Jacobi equation (EHJ) is obtained, and Eq. Q reads 



(20) 
(21) 

(22) 



The same equation is recovered directly when we introduce the transformation on the 
canonical momentas 



n 



as 



in Eq. ((7j), in the new coordinates x and y, this equation takes the form 



OS 
dx 



I) — . 



(23) 



(24) 



and choosing S = SxSy we obtain the following solutions 



-x/2 



Sy = /i = cte. 



(25) 



so we get. 



Equation 



S(x,y) = ±2v^e--/2. 
in the new variables become 



6 

VP 

X 

7f 



-x/2 



-x/2 



6 
A 

VP 



,3a- A,/, 



(26) 



(27) 
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The classical behaviour is considered solving the relations between fj27|) and Eqs. ()5al5bp . 
obtaining the relation for the the scale factor and the scalar field 

a = aoe^*, (28) 

the corresponding time behaviour (dr = Ndt), 

2 

a = a.QT^ , 

2 / A2 



having an inflationary scenario for the scale factor, with increasing power law when A < -\/2- 
This is the result from the classical approach (Chimento and Jakubi,1996). Unfortunately 
this potential has a serious drawback, once we fix the value of A to for inflation, there is no 
way to end it. 



III. SUPERSYMMETRIC FACTORIZATION SCHEME 

As already mentioned, the goal of this paper is to use the factorization approach of 
supersymmetric quantum mechanics to obtain the family of isoespectral potentials to V{(f)) = 
e'^"^ , an see if these potentials lead to new physics. 

If a quantum system is characterized by Hamiltonian H with eigenvalues z/, we may ask 
if there exists another Hamiltonian (i.e. different potential) which has the same spectrum. 
This is an interesting question because we may not be able to distinguish one system from the 
other by making them physically equivalent. With this in mind, we will begin by reviewing 
the factorization approach. 

Lets start with the hamiltonian that appears in the study of scalar fields interacting with 
gravitation Eq. ljlH IT2|l . written in generic form as 

- ^ + Vi(qi)Pi = EiPi, Pi = (X,Y), qi = (x,y). 



Ei 



2, ,2 ' 



6_ X 

.,2 C 



Vi = <^ . (30) 
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It is easy to show that the first order differential operators 

•At = -^ + Wi(q'), (31) 

A^ = ^ + Wi(q'), (32) 

factorize the hamiltonian, l.h.s.of Eq. ()H()|l (here Wi plays the role of a superpotential func- 
tion) as 

H+ - Ei = A^A; i = 1, 2. (33) 

The potential term Vi(q^) is related to the superpotential function Wi(q'^) via the Riccati 
equation 



dWi 

Vi(q')-Ei = Wf-^, i = l,2 (34) 



Making the transformation 



u' 



Wi = — -, i = l,2 (35) 

Ui 

where the / means ()34|) is transformed into the original hamiltonians applied to the 
functions Ui that correspond to the solutions ()14I15|) . this implies that once we have a 
solution to the original Schrodinger like equation (jHUj) the superpotential function can be 
constructed. In this factorization scheme, Vf is the partner superpotential of Vf, and can 
be calculated by performing the product 

Hr - Ei = ArAt, H-fi(q^) = Eifi(q'), (36) 

where fi is the wave function related to the hamiltonian Hj". Then, the isospectral potential 
to V+(q^) is 

Vr(q')-Ei = Wf + W:. (37) 
Using we get the functional relation between Vj" and V^, being 

Vr(q')=V+(q^) + 2W:. (38) 
However, this is not the most general solution as will be shown in the following section. 



A. General Solution 



The general solution to the Ricatti equation (jHEI) is found from (Cooper, Khare and 
Sukhatme, 1995; Mielnik, 1984) 

Vr (q') - Ei = Wf + W; = Wf + Wj, (39) 
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which by choosing 

Wi = Wi + -, (40) 
Yi 

leads to a Bernoulh equation for yi, 

y:-2Wiyi = l, (41) 

whose solution is 

yi(q^) = H^, (42) 

where li(q') = Jq nfdx and 7j is the factorization parameter (7^ plays the role of a " time-like 
parameter" in the evolution of the isospectral wavefunction). 
By using Eq. lj^ . we may write Eq. (pn|) as 

Wi(qO = Wi + (43) 
li + 7i 

and, the entire family of bosonic potentials can be constructed from 

V+(q\ 7i) - Ei = Wf (q', 7i) - W[{q\ 7i), (44) 



or 



V+(q\7i) = Vr-2W: 



VrW)-4^.2^, (45) 



finally 

is the isospectral solution of the Schrodinger equations (|14ll5p for the new family of po- 



Ui = g(7i)f^, (46) 
li + 7i 



tentials ()45|1 . with the condition on the function g(7i) = A/7i (71 + 1 ) , though in the limit 

7i ^ 00 g(7i) 7i and Ui Ui. (47) 

Considering the particular case for A < V2 (which correspond to inflation), we plot the 
solutions in the variable x = 6a — X(j) for the isospectral potential ()45|) (figures 1, 2 and 3) 
and the corresponding wavefunctions ()46|) (figures 4 and 5). The total WDW isospectral 
wave function has the following form 



[Ji,(±2v/^e-/2)+J_i,(±2v/^e^/2)] 
*iso(x, y; 71, 72) = 



aoe'^^y + aie-'^^y 



Ii + 7i I2 + 72 



^) 
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Original potential Isopotential with g=2 
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Isopotential with g=5 Isopotential with g=10 





-15-10 -5 5 -15-10 -5 5 

Isopotential with g=50 Isopotential with g-100 




-15-10 -5 5 




-15-10 -5 5 



FIG. 1: Plot is of the original potential and the iso-potentials. At large scale, the original potential 
and the iso-potentials do not differ much. 
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FIG. 2: At small scale the iso-potential have new structure for small values in the parameter 7, 
giving rise to a mechanism to end inflation that is not present in the original potential. 

The 7i parameters are included not for factorization reasons (these wave functions, in quan- 
tum cosmology are still non normalizable, except when a wave packet is constructed) , but as 
decoherence parameter embodying a sort of quantum cosmological dissipation (or damping) 
distance. The wave function is highly oscillatory, but when the value of 7j is increased. 
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Original potential 



^-™^0 -30 




-20 



-10 



FIG. 3: The corresponding changes in the structure for the isopotential, mainly for small values 
in the 7 parameter. 

the iso-spectral wave function tends to the wave function of the original potential (this is 
similar to the transitory effects in quantum mechanics, i.e. probability density in a potential 
barrier). Unfortunately this parameter ji can not be used to go from one physical state to 
another, so it is not a time parameter in the usual sense. In other words, this parameter 
seems to play the role of a "supersymmetric time- like parameter" that gives the evolution 
from a supersymmetric state to the original state, and the usual behavior for the wave 
function is reached (see figures 3 and 5). 

In the plots of the potentials, we see that the difference between the original potential 
and the isopotential are small at large scale (figure 1), but the corresponding wavefunctions 
present a drastically different behavior to a point that it vanishes for the limit 'ji ^ and 
recovers the original shape when ji 00. But in all cases, the roots of the wave functions 
remain the same. This can be attributed to the different structure of the iso-spectral poten- 
tial, as it appears in figures 2 and 3, in the sense that the amplitude of the wavefunctions 
corresponding to the iso-potentials in the regime a: < 0, is reduced dramatically. 

One expects that these effects should appear in the dynamics of our universe, in particular 
in the inflationary epoch. For inflation we need material with the unusual property of 
negative pressure (i. e. scalar flelds), as already mentioned, the potential V{(f)) = e~^'^ is 
one of the most studied. The standard approximation technique for analyzing inflation is 
the slow-roll approximation given by the slow roll parameter e = i(v^)^ << 1 in units 
of Mp = 1, in particular for the exponential potential e — << 1, this means that once 
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we fix A to have infiation, there is no mechanism to end it (we can not violate the slow 
roll condition). For the iso-spectral potentials we can see from the plot that we have an 
oscillatory function superposed to the original potential and this oscillatory behavior leads 
to mechanism that ends infiation. To see this, we start by writing the slow roll parameter 
in terms of the potential V{x), it takes the form e = ^(^^i^^F^ + where a = and 
b — -y^z^-i for the original potential = 1- So again, once we fix the value of A to 
have infiation there is no way to end it. For the iso-spectral potentials the approximate 
relation holds ■i7^~-i7^ where Viso{x) is the iso-spectral potential, and again we may 
write the slow roll parameter in a similar way e = ^a^( !;°°|^| + 6)^, where the constants a 
and h depend at most on A, now in contrast with the original potential, once we fix the 
value of A for inflation, the ratio can be as large as we want. We can see this from 

the plots of the iso-potentials. If 7j is small, the potential is oscillatory and has several 
roots. In particular near one of the roots Xi V-^^ ^ 0, but V{x) — > that is, 

in the neighborhood of Xi, ^^f^ is large, violating the slow- roll conditions, this gives the 
mechanism to end inflation!. Moreover, once inflation has ended the iso -potentials have 
several minima, where the scalar fleld can oscillate and start the reheating process. Finally 
there is a non zero probability that the scalar field may tunnel out of this minimum and 
give start another infiationary epoch. 



IV. CONCLUSIONS AND OUTLOOK 



In this paper we have applied the factorization method to infiationary cosmology, we 
started with a scalar field model coupled to and FRW model and an exponential potential 
V{(f)) = VqC^^'^. Wc have shown that in the variables x and y we can reproduce the known 
result, that for A < 2 and get power law infiation. This new set of variables is important 
to apply the factorization approach of SUSY-QM, from which we have obtained a complete 
family of iso-spectral potentials and their corresponding wave functions. We have shown in 
the plots, that we obtain the original potential and wave function in the limit 7^ — > oo. Out 
of this limit we have a large class of potentials that have the same quantum spectrum (this 
can clearly be seen in Fig. 5). Therefore all of these potentials describe the same quantum 
system. These iso-potentials have a very similar behavior for large scales, so at the classical 
level we won't see different dynamics. The drastic differences in the iso- wave functions 
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Original wavefunction Wavefunction with g=2 




Wavefunction with g=5 Wavefunction with g=10 




FIG. 4: The iso-wave function for different values for the 7 parameter. The behavior for small 
values of 7 can be seen. This different behavior can have effects in the quantum perturbations, 
and may modify structure formation in the large scale universe. 

Original wavefunction 




-3L 

FIG. 5: This plot show how the iso-wave functions tend to original wave function when the 
parameter 7 — ^ 00. 

and the iso-potentials, have an impact on the the dynamics of the universe. We have seen 
that this family of iso-potentials, have a natural mechanism to end inflation, and present 
the possibility of new inflationary epochs. We have found that the 7, parameter plays the 
role of a dccoherence parameter that embodies a type of quantum cosmological damping 
distance. In other words, it plays the role of a supersymmetric time-like parameter. This 
procedure can be applied to other models like, tachyon driven cosmology ( Garcia-Compean, 
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Garcia- Jimenez, Obregon and Ramirez, 2005), or to toy models of the landscape of string 
theory (Kobakhidze and Mersini-Houghton, 2004). These results may point to the existence 
of other potentials with the same eigenvalue spectrum, so that in principle the landscape 
may be falsifiable. These ideas are being explored and will be presented elsewhere. 
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